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Abstract 

Background, Aims and Scope. The traditional method of using 
Global Warming Potentials (GWPs) to assess the effects of cli¬ 
mate change in Life Cycle Assessment (LCA) does not account 
for indirect atmospheric effects like saturation effects and emis¬ 
sions of ozone precursors. The aim of this paper is to provide 
GWP values for LCA purposes of the most important climate 
related gases including indirect effects and assess whether these 
values are dependent of different background conditions and 
size of change in emissions fluxes. 

Methods. In order to reflect atmospheric interactions between 
direct and indirect greenhouse gases, an atmospheric climate 
and chemistry model needs to be used to determine GWP val¬ 
ues. Here, the IMAGE model is used to determine GWP values 
for a time horizon of 100 years. Different runs were performed 
to test the outcomes against various background emission sce¬ 
narios and various sizes in emission fluxes. 

Results and Discussion. The GWP values for C0 2 , N 2 0 and 
CH 4 depend on the chosen scenario and are lower under sce¬ 
narios with higher background emissions. The GWPs of Halons, 
PFCs and SF 6 are scenario and emission flux independent. The 
GWPs of HFCs increase with higher background emissions and 
higher additional emission fluxes, whereas GWPs of CFC are 
scenario independent, but change when other emission fluxes 
are applied. Finally, GWPs of HCFCs are higher in scenarios 
with higher background emission scenarios, and decrease when 
larger emission fluxes are applied. The GWP values calculated 
with IMAGE with direct effects only are comparable with IPCC 
values. Inclusion of the indirect atmospheric effects changes some 
of the values positively or negatively. For CFCs this results in a 
value up to 70% lower. For Halons this results in a value up to 
1000% lower, resulting in a strong negative GWP. Inclusion of 
indirect effects increases the GWP of CH 4 by 50%. Newly in¬ 
troduced greenhouse gases which have only indirect effects are 
given here a GWP. S0 2 has a negative GWP which depends on 
both the flux size and the chosen background scenario. CO and 


NMVOC have a positive flux-and-scenario dependent GWP. The 
GWP of NO x can, dependent on de chosen conditions, be posi¬ 
tive or negative. 

Conclusions. The GWP values of this paper are a first attempt 
to provide a consistent set of GWP values for all direct and 
indirect greenhouse gases, including differences in GWP values 
per background scenario and flux size. The inclusion of indirect 
effects in GWP values causes large differences in some impor¬ 
tant greenhouse gases, which should not be ignored in LCA- 
analyses. We suggest for LCA purposes to use GWPs including 
indirect effects for marginal change and for the most realistic 
background scenario. However, existing uncertainties in the in¬ 
direct effects of greenhouse gases demand for a better under¬ 
standing of the importance of these effects. 

Keywords: Climate change; global warming potentials (GWP); 
IMAGE; impact assessment; LCA; midpoint characterization 
factors 


Introduction 

It is the purpose of this commentary to discuss the sensitiv¬ 
ity of Global Warming Potentials (GWP) of the most im¬ 
portant climate-related gases towards the inclusion of indi¬ 
rect climate effects. We welcome critical comments to our 
proposal. 

Climate change is one of the impacts considered in Life Cycle 
Assessment (LCA) (Udo de Haes et al. 1999a,b, Jolliet et al. 
2004). The rapid increase of atmospheric greenhouse gases 
concentrations due to energy, industrial and land use emis¬ 
sions causes significant changes in the climate. Climate 
change affects ecosystems and possibly a wide range of hu¬ 
man health aspects (IPCC 2001a). The Global Warming 
Potential (GWP) is generally used to indicate the impact of 
climate change in LCA (Guinee 2001). This indicator gives 
information about the relative impact of greenhouse gases; 
it is a measure of the total increase in radiative forcing per 
unit of emission. GWP is a useful measure to compare and 
sum the effects of different greenhouse gases. Disadvantages 
include not accounting for geographical differences (Eickhout 
et al. 2003) and the arbitrary time horizon (Smith & Wigley 
2000a,b, Van Vuuren et al. 2006). In comparison with alter¬ 
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native methods it is however the most feasible method for 
political purposes (Fuglestvedt 2003). 

GWPs used for LCA (Guinee et al. 2001) are based on the 
list available from IPCC for direct greenhouse effect of the 
gases. Because of the great uncertainties about their magni¬ 
tudes, indirect effects were excluded in deriving the IPCC 
GWPs (IPCC 2001,2007). Indirect effects on the GWP arise 
from the formation and depletion of other atmospheric sub¬ 
stances due to an emission of a gas. Effects of the so-called 
indirect greenhouse gases CO, NO x , NMVOC and S0 2 were 
not considered in the IPCC-GWPs. These substances are no 
greenhouse gases themselves, but have an effect on the for¬ 
mation of tropospheric ozone, methane and sulphate aero¬ 
sols, which do influence the radiative forcing (IPCC 2001b). 
In addition, well-known greenhouse gases have indirect ef¬ 
fects as well. Chlorides such as chlorofluorcarbons (CFCs), 
and hydrochlorofluorcarbons (HCFCs) for instance, increase 
the radiative forcing, but are also related to the depletion of 
stratospheric ozone, causing a cooling of the atmosphere. 
This indirect effect can be even larger than the direct effects 
(Daniel et al. 1995, WMO 2003, Velders et al. 2005). For 
the purpose of LCA impact assessment, it may be worth¬ 
while to consider the indirect effects, in spite of the great 
uncertainties involved. 

Recent studies in the LCA field focused on the adaptation of 
existing risk assessment models to derive and improve 
characterisation factors for various impact categories, such 
as acidification (Seppala et al. 2006), photochemical ozone 
formation (Labouze et al. 2004), stratospheric ozone deple¬ 
tion (Hayashi et al. 2006), and toxicity (McKone & Hertwich 
2001). In this study we use a similar approach for the im¬ 
provement of characterisation factors for climate change. 
The Integrated Model for the Assessment of the Global En¬ 
vironment (IMAGE) simulates changes in climate variables 
and land use patterns as result of changes in demographic, 
technological and economic changes in society through 
changes in emissions of greenhouse gases (Alcamo 1994, 
IMAGE Team 2001). The model can consistently calculate 
the increase in radiative forcing as a result of emissions of 
well-known greenhouse gases (C0 2 , N 2 0, CH 4 , Hydrofluor- 
carbons (HFCs), Perfluorcarbons (PFCs) and SF 6 ), priority 
air pollutants (S0 2 , NO x , CO and non-methane volatile or¬ 


ganic compounds (NMVOC)) and ozone depleting pollut¬ 
ants (Halons and Chlorides). The IMAGE 2.2 model includes 
feedbacks from the ocean and the terrestrial environment to 
the atmospheric environment and vice versa (Alcamo 1994). 
The IMAGE model has been documented in detail (Alcamo 
1994, Alcamo et al. 1998, IMAGE Team 2001, MNP 2006) 
and discussed in scientific papers (e.g. Leemans & Eickhout 
2004, Strengers et al, 2004, Eickhout et al. 2006, Van Vuuren 
et al. 2007). Moreover, the IMAGE model has been used in 
several global assessments to address issues like greenhouse 
gas emissions (IPCC 2000), global environmental change 
(GE03 2002) and ecosystem goods and services (MA 2005). 
The latest scientific review of the IMAGE model concluded 
the general appreciation in policy circles has been built on 
the scientific credibility of the IMAGE model (Hordijk et al. 
2007). In the latest IPCC Fourth Assessment Report, results 
of the IMAGE model are also used for conclusions on cli¬ 
mate mitigation (IPCC 2007). 

This paper focuses on the atmospheric effect of greenhouse 
gases only. In IMAGE 2.2 the influence of four emission 
scenarios are assessed, reflecting different views on how the 
world may develop in the next 100 years. In the Special Re¬ 
port on Emission Scenarios (SRES) (IPCC, 2000) these so- 
called SRES scenarios are described. The goal of this paper is 
(i) to calculate GWPs of direct and indirect greenhouse gases 
including indirect atmospheric effects with the IMAGE 2.2 
atmospheric climate and chemistry model (Eickhout et al. 
2004) for use in Life Cycle Impact Assessment (LCIA) and (ii) 
test the linearity of these IMAGE-based midpoint character¬ 
ization factors against the chosen emission change applied 
and the chosen emission background scenario. 

1 Methodology 

1.1 IMAGE 

1.1.1 Model structure 

IMAGE 2.2 describes the climate change systems in all its 
complexity (Alcamo 1994). Although still a simplification, 
it includes many feedbacks to show the influence of most 
aspects of climate change. Fig. 1 shows a schematic over¬ 
view of the IMAGE model. 


r 



Fig-1 : Schematic overview of the interrelations of the three systems in the IMAGE model. Within the dotted line, the part of the model used for the study 
in this paper is selected 
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The model consists of three systems: 

• The Energy Industry System (EIS) computes all emissions 
of greenhouse and air polluting gases out of energy use 
and industrial activities with economic and population 
driving forces as main input of the system 

• The Terrestrial Environment System (TES) calculates food 
demand, land-use changes and changes in carbon uptake 
capacity, resulting in land use emissions of greenhouse and 
air polluting gases; it includes feedbacks of climatic change 
on the terrestrial environment, like C0 2 fertilization and 
the temperature effect on soil respiration 

• The Atmosphere Ocean System (AOS) computes radia¬ 
tive forcing, temperature change and sea-level rise due 
to changes in atmospheric concentrations of greenhouse 
gases as a result of the emission of gases out of the two 
previous systems 

1.1.2 Climate feedback 

Feedback from climate change also influences the system. 
The terrestrial uptake of carbon, calculated in the TES, is 
influenced by climate change. Furthermore, temperature 
change is used to calculate oceanic uptake of C0 2 . In this 
study we are interested in the effect of small additional emis¬ 
sions on radiative forcing. To get a clear view on how the 
individual substances behave atmospherically, the feedbacks 
of temperature on the terrestrial and oceanic system are left 
out of the calculation. Land use emissions and terrestrial 
carbon uptake are only considered exogenously (dependent 
per scenario) and oceanic carbon uptake is only dependent 
of the atmospheric C0 2 concentration determined by AOS. 
In Fig. 1, the parts of the model which are used are shown 
within the dotted line. More detail on how these processes 
are implemented in IMAGE is given in Eickhout et al. (2004). 


1.1.3 Calculation of concentration of greenhouse gases 

The calculation of the concentration of a greenhouse gas 
from the emissions is given by the equation: 


a\ghg 


A t 


~ CV GHG ' EghG ^ghg 


[ghg] 



where [GHG] is the atmospheric concentration of greenhouse 
gas GHG (in ppmv), E GHG the emission of the GHG in kg/yr, 
k G HGthe reac ti° n rate which is the reciprocal of the atmospheric 
lifetime (yr -1 ) and cv GHG a mass-to concentration conversion 
factor (ppmv/kg). The lifetimes of N 2 0, CFCs, Halons and 
PFCs and SF 6 are assumed to be constant. 

The lifetimes of the gases CH 4 , HCFCs and the HFCs are 
related to the OH abundance: 

^ GHG - k GHG+OH ■[OH] + k 

srat+bio (2) 


the loss rate to stratosphere and biosphere, which are as¬ 
sumed to be constant. A higher background concentration 
of OH radicals results in a lower lifetime of the chemical. 
This means that the chemical compound disappears faster 
out of the atmosphere. 

Calculation of concentration of C0 2 is more complex given 
the uptake of carbon by oceans and terrestrial biosphere. 


A[ C0 2 ] 
At 




The sinks (S CQ2 in kg/yr) are the C0 2 uptake by vegetation, 
calculated in the TES, and the C0 2 uptake by the oceans, 
calculated in AOS. Here, sink fluxes from TES are pre-cal- 
culated and used as input exogenously and are not part of 
the simulation itself. The values of lifetimes, mass-to con¬ 
centration factors and reaction rates are given in the sup¬ 
porting information. 


1.1.4 Calculation of concentration of tropospheric ozone and 
OH-radicals 

The indirect greenhouse gases CH 4 , NO x , CO and NMVOC 
form tropospheric ozone and OH-radicals. 


ln[OH] - P CH * In [CH 4 ] + P NO • E NO • P co • E co 

4 XX 

. P . f 

1 NMVOC ^ NMVOC 



n[trop.0 3 ] - Pch 4 ' [CH 4 ] + Pno x ’ ^ no x ’ ^co ’ ^co ( 5 ) 

’ ^NMVOC ’ E N MVOC 

P and R in these formulas are empirically derived sensitivity 
coefficients to determine the change in OH abundance and 
tropospheric ozone concentration (IPCC, 2001b). The val¬ 
ues are given in the supporting information. Mind that the 
values of P are, except for NO x , negative while the values of 
R are positive. This means that with increasing concentra¬ 
tions of CH 4 , CO and NMVOC the OH abundance decreases 
while the tropospheric ozone concentration increases. Natu¬ 
ral emissions from sources like natural vegetation, volca¬ 
noes and oceans are included, although these are kept con¬ 
stant over time. Tropospheric ozone is a greenhouse gas. 
OH-radicals change lifetimes and thereby the concentration 
of the gases CH 4 , HCFCs and the HFCs according to Eq. 2. 

1.1.5 Calculation of change in radiative forcing 

The radiative forcing due to halocarbons, tropospheric ozone 
and SF 6 is modelled as proportional to the calculated con¬ 
centration (ppmv). The proportionality constant is the em¬ 
pirically derived radiative efficiency coefficient re GHG (Wnr 
2 ppbv _1 ). The gas-specific re GHG values used are given in the 
supporting information. 


where k GHG+OH represents the oxidation of the GHG with 

OH, [OH] the concentration of OH-radicals and k strat+bio ^ ghg, t = re GHG ’ L GHG \ t ( 6 ) 
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The radiative forcing of C0 2 depends on pre-industrial con- 
centrations: 





\c02V 

lC0 2 ],o , 



The radiative forcing due to CH 4 and N 2 0 is besides the 
pre-industrial concentration of the substances themselves also 
dependent on the pre-industrial concentration of each other. 
The equation of the RF of both substances is given as: 

RFgcg =re GHG J[GHG\ t ~^[GHG\ i(< 

( 8 ) 

-f{[GHGl,\GHG\ t0 ,\GHG2] t0 ) 


GHG is the greenhouse gas CH 4 or N 2 0 and GHG2 is the 
other greenhouse gas; in the equation of CH 4 this is N 2 0 
and in the equation of N 2 0 this is CH 4 . The function in Eq. 
8 is given in the supporting information. Also the values of 
the pre-industrial concentrations of C0 2 , CH 4 and N 2 0 are 
given there. 

The negative radiative forcing due to ozone depletion of 
chlorinated chlorides and bromated halons in the strato¬ 
sphere is represented by: 


RF. 


aero -indirect ,t 


logfl + E S 0 2 ,,IE S02 _ nat ) 

l°g(l + E S q 2 1990 / AvO, 

—nat ) 


■RF 


aero-indirect ,1990 



The initial emission of S02 in 1990 assumed is 71.6 Tg S. 
The initial radiative forcing in 1990 due to sulphate aero¬ 
sols are -0.3 Wm-2 and -0.8 Wm-2 for direct and indirect 
effects respectively. The natural emission of S0 2 is assumed 
to be a constant amount of 34.4 Tg S. 


1.2 Global warming potentials 

IPCC (2001b) distinguishes two different measures for glo¬ 
bal warming: the Absolute Global Warming Potential 
(AGWP) and the Global Warming Potential (GWP). We de¬ 
fine the AGWP as the increase over time in radiative forcing 
per kilogram emission change of a substance over a given 
time horizon: 


AGWP(x) = 


ARF th 

E x 



The nominator represents the change in radiative forcing at 
time TH. E x is the yearly added emission flux. The AGWP is 
expressed in W nr 2 kg -1 yr. The Global Warming Potential 
is the AGWP of a substance related to the AGWP of a refer¬ 
ence gas, C0 2 (C0 2 -eq.): 


RF S trat.0 3 ,t ~ 0.001 


d a x{NCl r [GHG]J- 


+ d Bl .'Z(NBr r [GHG ] u ) 


(9) 


where GHG is the concentration of the concentration of the 
chlorinated or bromated GHG and NC^ and NBq are re¬ 
spectively the number of reactive chlorine or bromine at¬ 
oms in GHG i. d cl and d Br are the depletion coefficients of 
chlorine and bromine. These numbers are given in the sup¬ 
porting information. 

The radiative forcing due to stratospheric water vapour is 
positively affected by a higher concentration of CH 4 : 

RFstmt. H 2 o,t =0.05- re CH4 • {j[CH 4 \ t -^[CH 4 \ t0 ) (10) 

Another forcing effect is the direct radiative forcing of tro¬ 
pospheric sulphate aerosols (Myhre 1998). This is a result 
of the emission of the indirect greenhouse gas S0 2 . 


GWP(x) = 


AGWP ghg 

AGWP CQo 



The approach used by IPCC (2001b), is a so-called pulse-based 
method. This means that at time zero an emission pulse of a 
substance is released, which results in an instantaneous in¬ 
crease in the concentration in the atmosphere, whereafter the 
concentration gradually decreases over time due to reactions 
in the atmosphere. The increase in concentration corresponds 
with an increase in radiative forcing. The approach used in 
this paper is the flux-based method. In this case, the extra 
emission is added yearly, resulting in an increase of the radia¬ 
tive forcing towards a new steady state situation. The situa¬ 
tion is shown in Fig. 2. Fluxes and pulses in a linear system 
yield the same characterization factor (Heijungs 1995). 

From the year 2000 an extra emission flux has been added 
to the standard scenario up to the year 2100. The time hori¬ 
zon of 100 years was chosen because it is common in cli¬ 
mate change research, as it is the average atmospheric life¬ 
time of C0 2 (IPCC, 2001b). GWPs were calculated for C0 2 , 
N 2 0, CH 4 , NO x , S0 2 , CO, NMVOC, the halocarbons 
(CFCs, HCFCs, HFCs and Halons), the PFCs and SF 6 . 


RF 


'S0 2 ,t 


aero-direct,t 


'S0 2 ,1990 


■RF 


aero-direct,1990 



Sulphate aerosols cause also indirect effects. Indirect effects 
of tropospheric aerosols include the changing of the cloud 
albedo and amount due to cloud condensation and ice nu¬ 
clei. For this effect also the natural stratospheric emission of 
S 0 2 Es 02 -nat is considered: 



Fig. 2: Change in RF with the flux-based method 
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1.3 Simulations 

In the evaluation of GWPs with IMAGE several choices were 
required. To assess the influence of the percentage of mar¬ 
ginal change and the choice for a background emission sce¬ 
nario, the following simulations were performed. 

1.3.1 Marginal emission change 

We tested the linearity of the effect of change of emissions 
on the change in radiative forcing. The simulations were 
executed for different amounts of AE. We tested this linear¬ 
ity with respectively 1%, 10% and 100% extra emission of 
the total anthropogenic emissions in the year 2000. In a lin¬ 
ear situation the RF will increase linearly with the emission 
what will lead to a GWP value which is the same for all the 
three emission changes. If this is not the case the GWP val¬ 
ues for the three emission changes will be used to calculate 
the GWP that reflects a marginal emission change. This is 
done by making an exponential fit, shown in Fig. 3. The 
marginal change is the change closest to zero. 



Fig. 3: The GWP for different scenario fluxes can be exponentially re¬ 
lated. The GWP for the marginal emission change is the value for the 
smallest possible emission 


1.3.2 Background emission scenarios 

Marginal changes were made in the four different base sce¬ 
narios presented in IMAGE 2.2. The Intergovernmental Panel 
on Climate Change (IPCC) uses these four so-called SRES 
scenarios Al, A2, B1 and B2 (IPCC, 2000). IPCC developed 
these multiple scenarios to explore the uncertainties behind 
potential trends in global developments and GHG emissions. 
The scenarios differ in how global regions interrelate, how 
new technologies diffuse, how regional economic activities 
evolve, how protection of local and regional environments is 
implemented, and how demographic structure changes (IPCC 
2000). These different developments result in differences in 
projected emissions. These emissions are given in Fig. 4. 

2 Results and Discussion 

2.1 Direct GWPs 

Table 1 and 2 show direct GWPs of greenhouse gases. 'Di¬ 
rect' here means that only the change in radiative forcing 
due to changes in the concentration of the substance itself is 
taken into account. Indirect effects, such as the change in 
radiative forcing due to changes in tropospheric ozone and 
lifetime changes of CH 4 , HCFCs and HFCs and indirect forc¬ 
ing effects, are considered in the next result section. The 
reference value was the AGWP of 1 % C0 2 change in the Al 
scenario, being 9.06.10 -14 Wnr 2 kg _1 yr. 

2.2 Sensitivity of direct GWPs to the chosen emission change 

Table 2 shows the direct GWPs for the emission changes of 
1 %, 10% and 100% for the three main substances C0 2 , 
CH 4 and N 2 0. Larger emissions cause lower GWP values 
presenting a lower increase in radiative forcing per kilogram 
of added substance. Eq. 1 and 6, show that the change in 
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Fig. 4: The emission development in the 4 different scenarios. The emission development is given in relation to the year 2000 
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Table 1: Direct GWP values for the main greenhouse gases for different 
background scenarios and a 1% emission change relative to the total 
emission of the substance in the year 2000. Only direct effects of the 
substance on the radiative forcing are considered 


Substance 

Al 

A2 

B1 

B2 

C0 2 

1 

0.9 

1.5 

1.2 

n 2 o 

276 

268 

288 

281 

ch 4 

18 

15 

20 

18 

CFC11 

4,800 

4,800 

4,800 

4,800 

CFC12 

11,000 

11,000 

11,000 

11,000 

CFC113 

6,200 

6,200 

6,200 

6,200 

CFC114 

10,000 

10,000 

10,000 

10,000 

CFC115 

7,500 

7,500 

7,500 

7,500 

ecu 

1,800 

1,800 

1,800 

1,800 

CH3CCI 3 

2,100 

2,300 

2,100 

2,200 

HCFC22 

2,500 

2,700 

2,500 

2,600 

HCFC123 

150 

170 

150 

150 

HCFC124 

700 

900 a 

800 a 

800 

HCFC141b 

880 

960 

860 

910 

HCFC142b 

3,000 

3,300 

3,000 

3,100 

HCFC225ca 

700 

800 

700 

800 

Flalon1211 

1,300 

1,300 

1,300 

1,300 

Flalon1301 

7,200 

7,200 

7,200 

7,200 

CH 3 Br 

3 

3 

3 

3 

HFC23 

13,100 

13,200 

13,100 

13,100 

HFC32 

900 a 

900 

800 

900 

HFC43-10 

5,800 

5,600 

5,000 

5,200 

HFC125 

5,300 

5,200 

4,900 

5,000 

FI FC 134a 

1,900 

2,100 

1,900 

2,000 

FI FC 143a 

6,100 a 

5,600 

5,400 

5,500 

HFC227ea 

4,900 

5,300 a 

5,000 

5,200 

HFC245ca 

970 

1,080 

950 

1,000 

LL 

o 

6,100 

6,100 

6,100 

6,100 

c 2 f 6 

14,700 

14,700 

14,700 

14,700 

sf 6 

23,400 

23,400 

23,400 

23,400 


a HCFC124, HFC32, FIFC143 and FIFC227 have relatively low emis¬ 
sions. For small emission changes of 1% it is possible that due to 
rounding errors the GWP values differ slightly from values calculated 
with larger fluxes 


Table 2: The influence of different emission changes on the GWP of C0 2 , 
N 2 0 and CH 4 


Substance 

1% 

10% 

100% 

C0 2 

1 

0.98 

0.82 

ch 4 

28 

27 

25 

n 2 o 

276 

275 

265 


concentration and radiative forcing due to concentration of 
most substances is related to a constant value, while for C0 2 , 
CH 4 and N 2 0 (see Eq. 7 and 8) the relation between initial 
and actual concentration is involved. The lower GWP-val- 
ues by larger contributions can be declared by the fact that 
the contribution to global warming is less per molecule when 
there are more molecules present in the atmosphere. This 
saturation effect happens as the forcing gradually saturates 
into absorption bands while the strongest absorption bands 
have already been taken (Eickhout et al. 2004). This effect 
was also found by Caldeira & Kasting (1993). 


2.3 Sensitivity of direct GWPs for the chosen background 
scenario 

Table 1 shows that background emissions play a role in the 
effect of contributions of C0 2 , CH 4 , N 2 0, HCFCs, CEECCT 
and HFCs irrespective of indirect influences. The scenarios 
with the highest background emissions (see Fig. 4), A2 and 
Al, give the lowest GWP values for C0 2 , CH 4 and N 2 0. 
This is the same type of result as seen from the assessment 
of different emission changes: higher background emissions 
of these substances results in a lower GWP per added kilo¬ 
gram of the substance. For HFCs, CH 3 CC1 3 and HCFCs the 
scenarios with the high background concentrations result in 
higher values of the GWP. These gases do not have a con¬ 
stant lifetime but one that depends on the concentration of 
OH-radicals (see Eq. 2). The OH abundance decreases with 
higher emissions of CH 4 , NO x , CO and NMVOC (see Eq. 
5). This causes a higher lifetime for the HCFCs, CH3CC1 3 
CH 4 , and the HFCs, which results in a higher concentration 
growth and therefore a higher RF (IPCC, 2001b). CFCs, 
Halons and PFCs (including SF 6 ) are independent of back¬ 
ground scenarios. 

2.4 Indirect GWP effects 

Table 3 shows the GWPs of the greenhouse gases including 
indirect effects. The reference value was the same as the one 
used in Table 1 and 2 (9.06.10 -14 Wm -2 kg _1 yr). As can be seen 
in the table some substances have additional radiative forcing 
effects besides the direct effects shown in Table 1 and 2. Addi¬ 
tional effects are seen for CH 4 , CFCs, HCFCs and halons. 
The ozone precursors NO x , CO, NMVOC and S0 2 did not 
have direct effects, but they do show indirect effects. 

The values of the total GWPs are significantly lower for CFCs 
and the halons. This is a result of the depletion of strato¬ 
spheric ozone, which leads to a negative contribution to the 
radiative forcing resulting in a net lower GWP. As can be 
seen by Eq. 9 this effect caused by chlorine atoms causes an 
exponential effect in the change in radiative forcing. The bro¬ 
mine atoms do not have this exponential effect. The GWP- 
values for CFCs for small emissions are about 30% lower 
then the ones including only direct effects. With additional 
emissions up to 100% (see Table 3) the GWP values can even 
be 70% lower. In the case of halons the indirect cooling effect 
is about 10 times higher then the direct warming effect, result¬ 
ing in a negative GWP value. This results in a strong net nega¬ 
tive effect of these substances. The indirect effect for HCFCs 
is generally low, about 1 %. These outcomes were also found 
by Daniel et al (1995). Interestingly, the combination of ef¬ 
fects for HCFCs leads to lower GWP values in scenarios 
with low background emissions due to more OH abundance 
(see Section 2.3), but higher GWP values with higher emis¬ 
sion fluxes due to more stratospheric depletion. 

Substances like CH 4 , CO, NMVOC and NO x have multiple 
indirect effects. They include the formation of the green¬ 
house gas tropospheric ozone and the OH radical. The last 
effect affects the lifetime of CH 4 , HCFCs and HFCs. CH 4 , 
CO, NMVOC to cause positive contributions through both 
effects, while NO x has a positive effect through tropospheric 
ozone but a negative effect through OH abundance. This 
can result in either a net negative or a net positive GWP 
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Table 3: GWPs of with IMAGE 2.2 calculated greenhouse gases including direct and indirect effects (substances with only direct effects are left out of this table) 



Al 



A2 



B1 



B2 



Substance 

1% 

10% 

100% 

1% 

10% 

100% 

1% 

10% 

100% 

1% 

10% 

100% 

ch 4 

28 

27 

25 

21 

20 

19 

31 

31 

27 

26 

26 

24 

CO 

2.7 

2.7 

2.8 

2.8 

2.8 

2.8 

2.4 

2.4 

2.4 

2.6 

2.6 

2.6 

NOx 

-1.2 

-1.1 

-0.4 

-1.3 

-1.2 

-0.9 

2.5 

2.6 

3.1 

0.0 

0.1 

0.5 

NMVOC 

8.2 

8.2 

8.2 

8.2 

8.2 

8.2 

7.2 

7.3 

7.3 

7.9 

7.9 

7.9 

S0 2 

-94 

-90 

-71 

-59 

-58 

-52 

-118 

-112 

-82 

-101 

-97 

-74 

CFC11 

3,200 

3,000 

1500 

3,200 

3,000 

1,500 

3,200 

3,000 

1,500 

3,200 

3,000 

1,500 

CFC12 

6,000 

5,700 

2,700 

6,000 

5,700 

2,700 

6,000 

5,700 

2,700 

6,000 

5,700 

2,700 

CFC113 

4,700 

4,600 

3,600 

4,700 

4,600 

3,600 

4,700 

4,600 

3,600 

4,700 

4,600 

3,600 

CFC114 

9,800 

9,700 

9,700 

9,800 

9,700 

9,700 

9,800 

9,700 

9,700 

9,800 

9,700 

9,700 

CFC115 

7,400 

7,300 

7,300 

7,400 

7,300 

7,300 

7,400 

7,300 

7,300 

7,400 

7,300 

7,300 

ecu 

1,300 

1,200 

640 

1,300 

1,200 

640 

1,300 

1,200 

640 

1,300 

1,200 

640 

CH 3 CCI 3 

2,000 

1,900 

1,000 

2,300 

2,100 

1,000 

2,000 

1,800 

1,000 

2,100 

1,900 

1,000 

HCFC22 

2,500 

2,400 

2,000 

2,700 

2,600 

2,200 

2,400 

2,400 

2,000 

2,500 

2,500 

2,100 

HCFC123 

150 

150 

150 

170 

170 

170 

150 

150 

150 

150 

150 

150 

HCFC124 

700 

700 

700 

900 a 

800 

800 

800 a 

700 

700 

800 

800 

800 

HCFC141 b 

870 

850 

740 

960 

930 

800 

860 

840 

730 

900 

880 

760 

HCFC142 b 

3,000 

3,000 

2,900 

3,300 

3,200 

3,100 

3,000 

2,900 

2,800 

3,100 

3,000 

2,900 

HCFC225 

700 

700 

700 

800 

800 

800 

700 

700 

700 

700 

700 

700 

Flalon1211 

-12,000 

-12,000 

-12,000 

-12,000 

-12,000 

-12,000 

-12,000 

-12,000 

-12,000 

-12,000 

-12,000 

-12,000 

Flalon1301 

-62,000 

-62,000 

-62,000 

-62,000 

-62,000 

-62,000 

-62,000 

-62,000 

-62,000 

-62,000 

-62,000 

-62,000 

CH 3 Br 

-824 

-824 

-824 

-824 

-824 

-824 

-824 

-824 

-824 

-824 

-824 

-824 


a HCFC124, HFC143 and FIFC227 have relatively low emissions. For small emission changes of 1% it is possible that due to rounding errors the 
GWP values slightly differ from values calculated with larger fluxes. 


depending on the chosen background conditions. CH 4 causes, 
besides these effects, also a forcing by water vapour (see 
Eq. 10). S0 2 causes the formation of sulphate aerosols. These 
particles have direct and indirect forcing effects (see Eq. 11 
and 12). Both effects cool the atmosphere. The greenhouse 
gases C0 2 N 2 0, the HFCs and the PFCs are found to have 
no indirect contributions and are therefore left out in Table 3. 

2.5 Direct and total GWPs in relation to previous research 

Table 4 shows the GWPs of the gases of scenario A1 with 
their direct and their total effects from Table 1 and Table 3 
for marginal emission changes. These values are compared 
with the values of IPCC (2001b). The table shows the GWPs 
for marginal contribution. For most greenhouse gases the 
GWP for 1, 10 and 100% emission change are approxi¬ 
mately equal. As stated in section 2.2 C0 2 , CH 4 , N 2 0, 
NMVOC, CO, NO x and the chlorides are dependent on the 
chosen emission change. Generally the GWP for marginal 
change is approximately the same as the GWP for 1 %. The 
marginal-change GWP of CFC11 and CFC12 is however 
smaller then the one for 1%. For these substances the expo¬ 
nential fit is used as given in section 1.3.1. IPCC uses a ref¬ 
erence value based on the Bern carbon cycle model assum¬ 
ing a constant future concentration of C0 2 . This value is 
8.99 10- 14 W nr 2 kg (IPCC 1994, 2001b). 

The direct GWPs calculated in this study are close to the 
values listed in the IPCC report (2001b). Noticeable differ¬ 
ences can be seen for CH 4 , chlorides, halons and HFCs. In 


the calculation of CH 4 IPCC did account for the indirect 
effects formation of stratospheric water vapour and tropo¬ 
spheric ozone. The effect of changing lifetimes was, how¬ 
ever, not included. An increase in lifetime due to a lower 
OH concentration (see Eq. 1, 2 and 4) affects the GWP also 
positively. This makes that the IPCC-value for CH 4 lies in 
between the value of direct and total GWP. The calculated 
GWPs for NO x CO, NMVOC and S0 2 were not listed be¬ 
fore in the IPCC report. However, Fuglestvedt et al (1996) 
investigated the indirect effects of CH 4 , CO and NO x through 
changes in tropospheric ozone and CH 4 resulting in esti¬ 
mated GWP values of 30 for CH 4 , 3 for CO and 'slightly 
negative' for NO x from surface sources. These results com¬ 
pare well with the ones found in this study. 

The results of the GWP values calculated with IMAGE in¬ 
cluding direct effects only are comparable to GWPs of IPCC 
(2001b). The total GWP values found in this study, which 
includes indirect effects of chlorides and halons, are much lower 
and compare well with results obtained by Daniel et al. (1995), 
WMO (2003) and Velders et al. (2005), who did similar stud¬ 
ies. However, since these studies used different emission sce¬ 
narios it is difficult to directly compare the results. 

The GWP values for HCFCs and HFCs differ the most from 
the values of IPCC. This can be seen as the effect of the way 
OH radical reaction rates are implemented in the GWP cal¬ 
culations. In IMAGE the lifetime of these GHGs explicitly 
depends on the presence of OH-radicals. The concentration 
of this radical is variable and is influenced by the emissions 
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Table 4: GWPs of with IMAGE 2.2 calculated greenhouse gases includ¬ 
ing indirect effects. In the last column the values of the GWPs as pub¬ 
lished in the report of IPCC (2001b) are shown 



GWP 100 direct 

GWP 100 total 

IPCC, 2001 

o 

o 

ro 

1 

1 

1 

n 2 o 

276 

276 

296 

ch 4 

18 

28 

23 

CO 

0 

3 


NOx 

0 

-1 


NMVOC 

0 

8 


S0 2 

0 

-94 


CFC-11 

4,800 

3,300 

4,600 

CFC-12 

11,000 

6,100 

10,600 

CFC-113 

6,200 

4,700 

6,000 

CFC-114 

10,000 

9,800 

9,800 

CFC-115 

7,500 

7,400 

7,200 

CCI 4 

1,800 

1,300 

1,800 

MCF (CH 3 CCI 3 ) 

2,100 

2,000 

140 

HCFC-22 

2,500 

2,500 

1700 

HCFC-123 

150 

150 

120 

HCFC-124 

700 

700 

620 

HCFC-141b 

880 

870 

700 

HCFC-142b 

3,000 

3,000 

2,400 

HCFC-225 

700 

700 

180 

Halon-1211 

1,300 

- 12,000 

1,300 

Halon-1301 

7,200 

-62,000 

6,900 

CH 3 Br 

3 

-824 

5 

HFC-23 

13,100 

13,100 

12,000 

HFC-32 

1,000 a 

1,000 a 

550 

HFC-43-10 

5,800 

5,800 

1,500 

HFC-125 

5,300 

5,300 

3,400 

HFC-134a 

1,900 

1,900 

1,300 

HFC-143a 

6,200 a 

6,200 a 

4,300 

HFC-227ea 

4,900 

4,900 

3,500 

HFC-245ca 

970 

970 

640 

'n|- 

LL 

0 

6,100 

6,100 

5,700 

c 2 f 6 

14,700 

14,700 

11,900 

sf 6 

23,400 

23,400 

22,200 


a For HFC32 and HFC143 and HFC227 the GWP of 10% is taken 
because the values for 1 % slightly differ due to rounding errors. 


of CH 4 , CO, NO x and NMVOC, while in the calculation of 
IPCC (2001b) constant atmospheric OH concentration was 
assumed. With the assumption that the emissions of CH 4 , 
CO, NO x and NMVOC will increase in the next century 
(see Fig. 4) the calculated lifetimes of IMAGE (see Eq. 2) 
will be higher resulting in a higher GWP. 

2.6 Normalisation 

In order to show the relevancy of our results, we caclulated 
sub-specific and total normalisation factors in C0 2 equiva¬ 
lents for global warming, base on greenhouse gas emission 


Table 5: Total C0 2 equivalents (kg) for the year 2000 calculated with GWPs 
from Table 3 for direct and total effects 



Emission 2000 
(kg) 

C0 2 eq. direct 
(kg) 

C0 2 eq. total 
(kg) 

C0 2 

2.85 10 13 

2.85 10 13 

2.85 10 13 

N 2 0 

1.15 10 10 

3.16 10 12 

3.16 10 12 

ch 4 

2.99 10 11 

5.38 10 12 

8.37 10 12 

CO 

1.04 10 12 

0 

3.13 10 12 

NO x 

1.16 10 11 

0 

-1.16 10 11 

NMVOC 

1.67 10 11 

0 

1.33 10 12 

S0 2 

1.30 10 11 

0 

- 1.22 10 13 

CFC-11 

4.06 10 ° 7 

1.95 10 11 

1.34 10 11 

CFC-12 

1.01 10° 8 

1.11 10 12 

6.14 10 11 

CFC-113 

3.86 10 ° 6 

2.39 10 10 

1.81 10 10 

CFC-114 

2.07 10 ° 6 

2.07 10 10 

2.03 10 10 

CFC-15 

8.73 10 ° 5 

6.54 10 ° 9 

6.46 10 ° 9 

CCI 4 

4.17 10 ° 5 

7.51 10 ° 8 

5.42 10 ° 8 

MCF (CH 3 CCI 3 ) 

3.57 10 ° 5 

7.49 10 ° 8 

7.13 10 ° 8 

HCFC-22 

3.00 10 ° 8 

7.51 10 11 

7.51 10 11 

HCFC-124 

3.93 10 ° 6 

2.75 10 ° 9 

2.75 10 ° 9 

HCFC-141b 

1.66 10° 8 

1.46 10 11 

1.45 10 11 

HCFC-142b 

5.09 1 0 07 

1 . 5310 " 

1 . 5310 11 

Halon-1211 

4.82 10° 6 

6.26 10° 9 

-5.78 10 10 

Halon-1301 

9.26 10° 5 

6.67 10° 9 

-5.74 10 10 

HFC-125 

7.40 10° 6 

3.92 10 10 

3.92 10 10 

HFC-134a 

1.30 10° 8 

2.47 10 11 

2.47 10 11 

HFC-143a 

5.40 10° 6 

3 . 3510 10 

3 . 3510 10 

sf 6 

5.22 10° 6 

1.22 10 11 

1.22 10 11 

Total 


3 . 9910 13 

3.43 10 13 


Emissions are taken from Wegener Sleeswijk et al. (2007). GWPs for 
marginal emission changes of scenario Al were used to obtain these data. 


data for the world in year 2000. Emission data were taken 
from Wegener Sleeswijk et al. (2007). The results are given 
in Table 5. 

Assessing these results we see that the influence of indirect 
greenhouse gases, CO, NO x , NMVOC and S0 2 is great. 
The impact of CH 4 becomes greater and Halons appear as a 
cooling factor, when indirect effects are taken into account. 
The latter will surely decrease in the coming century be¬ 
cause of the expected decrease of these substances due to 
successful environmental policy (SRES 2000). 

2.7 Uncertainties 

Despite of the importance of indirect effects of greenhouse 
gases on the climate system, new accurate values of GWP 
cannot be determined yet due to uncertainties in these ef¬ 
fects. Fig. 5 shows the radiative forcings as IPCC reports 
them in its latest Assessment Report (IPCC 2007). Values 
of the well-mixed greenhouse gases are considered to be of 
high confidence on the level of scientific understanding 
(LOSU, IPCC 2007), whereas the indirect effects of aero¬ 
sols are considered of low confidence (IPCC 2007). Also 
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□ IPCC (20D7) 


□ IMAGE 2.2 
(Eickhout 
el al..2004) 


Fig. 5: Global and annual mean radiative forcing from 1750 to present 
according IPCC (2007) and according IMAGE 2.2 as presented in Eickhout 
et al. (2004). Only values that are reported in Eickhout et al. (2004) are 
taken from IPCC. IMAGE results are from 1765 to 1995; IPCC from pre¬ 
industrial to 2005 


the broad uncertainty range that is attached to the indirect 
effect of aerosols clearly shows the current level of under¬ 
standing. Forcings of aerosols are better understood than 
at the time of the Third Assessment Report (IPCC 2001) 
due to improved in situ, satellite and ground-based mea¬ 
surements and more comprehensive modelling, but remains 
the dominant uncertainty in radiative forcing (IPCC 2007). 
Therefore, this understanding needs to be improved further, 
before definite values of GWPs including indirect effects can 
be applied in LCA studies. 


3 Conclusions 

Our GWP values are generally equal with the values given 
in the IPCC report, but do in addition include indirect at¬ 
mospheric effects such as the effects of changes in the ozone 
concentration and variable lifetimes. The inclusion of indirect 
effects in GWP values returns large differences in some impor¬ 
tant greenhouse gases and, therefore, is a valuable contribu¬ 
tion to LCA-analyses. The simulations showed that the GWP 
values depend on the amount of emission change and the cho¬ 
sen background emission scenario. This can cause variations 
in the net GWP value up to a factor 2 for the amount of emis¬ 
sion change and up to a factor 1.5 for background scenario. 
This variation is a point of interest for the use of GWPs. For 
LCA purposes where small emissions changes are considered 
the GWP for marginal contributions can be used. The user of 
the GWP can chose a background scenario, for instance Al 
because it is most commonly used, but needs to be aware 
that the variation between different background scenarios 
is up to a factor 1.5, especially for C0 2 . 
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Table SI: Lifetimes and Mass-to-concentration factors for Greenhouse 
gases with constant lifetimes 


Greenhouse gases 
with constant lifetimes 

Mass to 
concentration 
factor cv (ppb/Tg) 

Lifetime (yr); 1/k 

C0 2 

0.4688 a 

b 

n 2 o 

0.21 

120 

CFC-11 

0.04314 

45 

CFC-12 

0.04901 

100 

CFC-113 

0.03163 

85 

CFC-114 

0.03467 

300 

CFC-115 

0.03836 

1,700 

ecu 

0.03853 

35 

Flalon-1211 

0.03584 

11 

Halon-1301 

0.03984 

65 

CH 3 Br 

0.0350877 

0.7 

LL 

o 

0.068 

50,000 

C2F 6 

0.0508 

10,000 

sf 6 

0.041 

3,200 


a The unit of the Mass concentration factor of C0 2 is ppmv/kg C. 

b The lifetime of C0 2 is dependent on variable fluxes of C0 2 into the 
ocean and the biosphere. 
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Table S2: Reaction rates and Mass-to-concentration factors for 
Greenhouse gases with constant lifetimes 


Greenhouse gas 
with non-constant 
lifetimes 

Mass to 
concentration 
factor cv 
(ppb/Tg) 

Reaction rate 
kGHG+ 0 H of the 
oxidation by OH 
(cm 3 yr _1 ) 

Transport 
losses to 
stratosphere 
and biosphere 
(yr 1 ) 

ch 4 

0.00040 

1.274-10 -7 

18 

MCF (CH3CCI3) 

0.04689 

1.781-10 -7 

2,100 

HCFC-22 

0.07234 

8.482-10 -8 

2,500 

HCFC-123 

0.0409 

8.150-1CT 7 

150 

HCFC-124 

0.04583 

1.985-10 -7 

700 

HCFC-141b 

0.05349 

1.144.1CT 7 

880 

HCFC-142b 

0.06224 

5.773-10 -8 

3,000 

HCFC-225 

0.03082 

1.683-10 -7 

700 

HFC-23 

0.086 

4.640-10 -9 

13,100 

HFC-32 

0.116 

1.930-10 -7 

900 

HFC-43-10 

0.07442 

9.236-10“ 8 

5,800 

HFC-125 

0.05211 

3.552-10" 8 

5,300 

HFC-134a 

0.07442 

9.236-10“ 8 

1,900 

HFC-143a 

0.07142 

2.110-10 -8 

6,100 

HFC-227ea 

0.035 

2.830-10“ 8 

4,900 

HFC-245ca 

0.0448 

1.660-1CT 7 

970 
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Table S3: Sensitivity coeffiecients for the formation of OH radicals and 
tropospheric ozone 

Sensitivity P in Sensitivity R 

relation to the in relation to 

OH radical tropospheric 

ozone 

-0.32 6.7 

0.0042 0.17 

-0.000105 0.0014 

-0.000315 0.0042 

Table S5: Radiative efficiency coefficients for C0 2 , CH 4 and N 2 0 and 
their initial pre-industrial concentrations in 1765 

Radiative efficiency Pre-industrial 

coefficient re concentration 

5.325 (Wm- 2 ppmv _1 ) 278 ppmv 

0.12 (Wnrf 2 ppbv _1 ) 700 ppbv 

CH 4 0.036 (Wnrf 2 ppbv -1 ) 270 ppbv 

The calculation of radiative forcing due to methane and 
nitrous oxide (equation 8 in the main article) is given by 
the equations I and II. The values for a and b are 2.01 10~ 5 
and 5.31 10~ 15 respectively. 

RFch, = re CHt - VM , 0 - 

0A1 ln[l + a([CH 4 \[N 2 o]J J5 +b[CH 4 \([CH^[N 2 0\J- 52 ]- \ (1) 

0.47 4 + a([CH 4 ] (0 [N 2 o ] t0 )° 75 + b[CH 4 ] (0 ([CH 4 \ 0 [n 2 o\ 0 f 52 ]/ 




Green house gas 

C0 2 

N 2 0 


(Indirect) 
Greenhouse gas 

CH 4 

CO 

NO x 

NMVOC 


Table S4: Radiative efficiency coefficients for halocarbons, SF 6 C0 2 , CH 4 
and N 2 0 and the number of chlorine and bromine atoms per compound 

Green house gas 

Trop. 0 3 

CFC-11 _ 

CFC-12 
CFC-113 
CFC-114 
CFC-115 

CCI 4 _ 

MCF (CH3CCI3) 

HCFC-22 
HCFC-123 
HCFC-124 
HCFC-141b 
HCFC-142b 
HCFC-225 



a The unit of radiative efficiency of tropospheric ozone is V\T 2 DU _1 
(Dobson unit) Dobson units are used to express the concentration 
of tropospheric ozone. 


Radiative efficiency Number of chlorine 


coefficient 

(Wm" 2 ppbv _1 ) 


0.042 


or bromine atoms 
per compound 
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Expanded Damage Function of Stratospheric Ozone Depletion to Cover Major Endpoints 
Regarding Life Cycle Impact Assessment 

Kentaro Hayashi, Ai Nakagawa, Norihiro Itsubo and Atsushi Inaba 



Background. Stratospheric ozone depletion is one of the important 
environmental issues for LCIA. The National LCA Project of Japan 
has developed a framework of LCIA since 1998, which tackles the 
issue employing an endpoint approach. Although the basic compo¬ 
nents were available in 2000, it was required that the target end¬ 
points should be expanded in particular. 

Objective. This study aimed at expanding the scope of damage func¬ 
tion of ozone depletion in the LCIA framework of LIME. Damage func¬ 
tion gives potential and quantitative damage for each endpoint per 
unit emission of ODS. 

Methods. Marginal damage due to the unit emission of ODS was 
calculated for 13 substances for which quantitative information was 
available as follows: (1) the increase of UVB radiation at the earth’s 
surface per unit emission of ODS was estimated, (2) the increase of 
potential damage per unit increase of UVB radiation was estimated, 
(3) the increase of potential damage per unit emission of ODS was 


determined by connecting the two relationships, and (4) correcting by 
the atmospheric lifetime of ODS, so that the damage function was 
then obtained. For other ODSs regulated by the Montreal Protocol, 
their damage functions were estimated by multiplying the ratio of ODP 
compared to the corresponding reference substance by the damage 
function of this reference substance. 

Results and Discussion. The damage function of ozone depletion 
included the following endpoints: skin cancer and cataract for human 
health, crop production and timber production for social assets, and 
terrestrial NPP and aquatic NPP for primary productivity. And dam¬ 
age factors for each safeguard subject were also obtained. 
Conclusion. The damage function of ozone depletion could cover all 
ODSs regulated by the Montreal Protocol and also cover important 
endpoints. Uncertainty of damage function is also an important point 
to be elucidated. Preliminary studies of uncertainty analysis have be¬ 
gun for the damage function of ozone depletion. However, further 
analysis is required to comprehensively evaluate the uncertainty of 
the damage function. 
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